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DETECTING  DEMAND  CHANGES 


1.  Introduction 

One  of  the  foremost  problems  facing  the  Inventory  manager  is  that  of 
anticipating  changes  in  demand  patterns  and  adjusting  the  inventory  level 
accordingly.  If  future  demands  vere  exactly  known  then  of  course  there  is 
ho  problem,  f  or  •  then  the  inventory  level  can  be  adjusted  in  anticipation 
of  demand  changes.  The  realistic  case  arises  when  demand  is  probabilistic 
over  time,  that  is,  a  stochastic  process.  Even  when  one  is  willing  to 
assume  the  mean  value  function  of  such  a  process  as  known,  it  is  not 
possible  to  predict  the  exact  demands  and  the  best  one  can  do  is  provide 
reasonable  estimates  of  future  demands.  Then  it  is  quite  appropriate  to 
adjust  the  inventory  level  in  accordance  with  those  estimates  in  such  a 
way  as  to  effect  rapid  response  to  changes  that  may  occur  as  time  progresses. 

If  the  inventory  level  at  any  particular  time  is  established  as  a 
result  of  analyzing  the  history  of  demands  up  to  that  time,  then  it  is 
quite  natural  to  view  an  inventory  system  as  a  servomechanism  or  "black 
box"  in  which  the  input  is  demand  and  the  output  is  inventory.  Such  an 
approach  to  the  study  of  inventory  systems  appears  to  have  begun  with  Simon 
[4]  and  Vassian  [5]  in  the  early  1950's,  was  revised  somewhat  briefly  by 
Howard  [2]  in  1963  and  examined  more  recently  by  Reilly  [3]  and  the  writer. 
Except  for  a  few  isolated  articles  in  the  literature,  no  extensive  appli¬ 
cation  of  the  numerous  tools  developed  in  servo  theory  to  inventory  systems 
seems  to  exist.  The  purpose  of  this  report  is  to  direct  attention  to 
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these  techniques  as  a  means  of  analyzing  and  designing  inventory  control 
systems  . 

No  attempt  vill  be  made  in  this  preliminary  report  to  obtain  numerical 
solutions  to  a  particular  problem.  Rather,  attention  vill  be  devoted  to 
structuring  a  model  vithin  vhich  it  is  anticipated  that  several  special 
problems  can  be  solved.  In  particular,  the  basic  problem  of  designing  an 
inventory  system  that  vill  be  sensitive  to  demand  changes  vill  be  borne 
in  mind  at  all  times.  Of  special  concern  is  that  of  an  item  vhose  demand 
pattern  changes  from  that  of  a  slow  mover  to  a  moderate  or  fast  moving  item. 
Hopefully,  there  vill  be  an  opportunity  in  the  future  to  process  some  real 
demand  data  from  a  Naval  supply  system  to  test  the  model  and  possibly 
suggest  changes  in  either  the  model  or  the  supply  system  as  a  result.  In 
this  early  stage,  relatively  single  examples  vill  be  examined  for  further 
insight  and  for  illustrative  purposes. 

2.  Description  of  the  Model 

The  inventory  model  ve  wish  to  adopt  in  this  paper  is  a  system  vith 
periodic  reviev  and  fixed  order  points.  Such  a  system  is  videly  used  in 
practical  inventory  situations  vhere  periodic  checks  of  stock  are  far  more 
economical  than  continuous  reviev.  It  is  one  that  has  actually  been  used 
in  the  past  by  the  military  for  spare  parts.  For  the  sake  of  convenience, 
ve  assume  that  order  points  are  fixed  and  equally  spaced.  Without  loss  of 
generality,  let  the  order  points  be  designated  by  the  positive  integers  on 
the  time  scale.  We  suppose  a  physically  realizable  system  so  that  only 
the  non-negative  time  axis  is  appropriate. 
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fl. 

To  be  more  specific,  ve  assume  that  the  t  demand  period  is  measured 
by  real  time  betveen  t-1  and  t  for  t  =  1,2,3, ...  •  At  the  end  of  a 
demand  period,  and  only  then,  a  replenishment  order  is  placed  against 
future  demands.  We  suppose  a  fixed  lag  time  T  in  delivery  vhere  T  is 
measured  in  demand  periods  so  that  an  order  placed  at  time  t  is  delivered 
at  time  t  +  T  and  is  hence  usable  for  fulfilling  demands  during  the 
(t+T+l)st  demand  period.  We  further  suppose  that  any  oversupply  at  the 
end  of  a  given  period  may  be  returned  to  vendor  so  that  a  replenishment  order 
may  be  negative. 

As  to  the  assunptions  about  the  inventory  level,  ve  suppose  that  the 
system  begins  vith  a  safety  level  or  initial  inventory  I  at  time  zero. 
Since  our  primary  concern  vill  be  vhether  or  not  the  inventory  level  in 
future  reviev  periods  is  greater  than  or  less  than  the  safety  level  ve  shall 
henceforth,  for  the  sake  of  simplicity,  assume  that  the  inventory  level,  1^, 

i.u 

at  the  end  of  the  t  demand  period  is  the  difference  betveen  the  actual 
inventory  and  IQ  at  that  time.  Moreover,  if  demand  should  exceed  on  hand 
inventory  ve  allow  back  orders  to  be  placed  against  future  deliveries.  Thus, 
negative  values  of  1^  are  permissible  in  the  analysis  to  follov. 

As  to  the  demand  itself,  the  above  conditions  make  it  natural  to  assume 

that  it  can  be  described  by  a  discrete  parameter  stochastic  process.  We 

th 

will  denote  the  total  demand  during  the  t  period  by  X^.,  t  =  1,2,3, ••• 
and  ve  let  XQ  =  0  for  the  sake  of  completeness.  It  vill  be  helpful  in 
the  analysis  if  ve  suppose  that  X^.  =  m(t)  +  vhere  m(t)  is  the 

deterministic  mean  value  function  of  the  process  {X^.}  and  is  a  random 
variable  vith  Efe^]  =  0  for  all  t.  The  process  {e^}  is  sometimes 
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referred  to  as  the  noise  component. 

The  mathematical  model  for  the  inventory  system  that  we  have  described 
can  now  be  nicely  summarized  by  the  recursive  relation, 

(2-1)  It  =  It_1  +  -  Xt,  t  =  1,2,3,... 

Such  an  inventory  system  has  the  feature  that  it  is  automatically  controlled 
once  we  decide  on  the  form  of  the  reorder  rule  0^.  We  turn  our  attention 
to  this  matter  in  the  next  section. 

3.  Determination  of  the  Reorder  Rule 

There  are  of  course  infinitely  many  ways  to  specify  a  reorder  rule 
consistent  with  our  model.  To  choose  one  as  best  requires  some  criterion 
for  choice.  An  obvious  criterion  in  the  present  case  is  to  require  an  order 
rule  which  will  minimize  the  inventory  variations  about  the  safety  level. 

It  is  not  surprising  that  a  solution  among  all  possible  choices  of  0  cannot 
exist  and  we  must  restrict  our  order  rules  to  lie  in  some  special  class  of 
functions  if  we  hope  to  find  a  solution.  A  natural  restriction  is  to 
consider  only  those  rules  which  depend  upon  past  demands  and  inventory 
levels.  To  this  end,  we  restrict  our  attention  to  only  those  order  quantities 
that  can  be  expressed  as  a  linear  combination  of  past  demands  and  inventories. 
However,  as  time  progresses,  we  would  certainly  want  to  weight  the  early 
part  of  our  history  with  smaller  and  smaller  values.  This  suggests  that 
we  express  a  general  reorder  rule  as  follows, 
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(3-i)  *t  -  J  Vt-j  ♦  J  Vt-y  t  =  1's'3'--' 

The  problem  Is  thus  reduced  to  that  of  determining  the  sequences 
^^0,^1,^2#  •  •  •  }  ^0*^1*  v-)-  Observe  that  the  same  weights  are 

always  given  to  the  most  recent  observations  while  the  weights  given  the 
early  observations  vary  with  time.  Of  course  Xq  and  Iq  are  both 
zero  but  it  does  no  harm  to  write  them  in  the  expression  for  9 ^  and  makes 
for  a  certain  degree  of  completeness  in  the  sequel. 

Even  with  the  additional  restrictions  it  is  not  possible  in  general  to 

find  a  unique  order  rule.  Vassian  [5]  has  given  a  partial  solution  to  the 

problem  in  the  sense  of  providing  a  rule  which  depends  upon  one*s  ability 

to  forecast,  at  each  time  t,  the  total  demand  over  the  following  T+l 

periods,  i.e.,  demands  during  periods  t+l  through  t+T+1.  Indeed,  this 

is  about  as  close  to  a  solution  as  one  would  expect  to  arrive,  Vassian* s 

claim  for  a  solution  notwithstanding,  because  of  the  random  nature  of 

demand.  At  any  time  t,  past  demands  are  of  course  known  but  future  demands, 

hence  inventories,  are  random.  Consequently,  the  concept  of  minimizing 

inventory  should  be  replaced  with  that  of  minimizing  expected  inventory  and 

we  require  that  at  least  lim  E[I.  ]  =  0  for  any  specified  order  rule. 

t  -+  oo 

It  is  then  possible  to  find  a  family  of  order  rules  which  will  satisfy 
equations  (2-l)  and  (3-l)  subject  to  this  condition.  The  details  of  the 
solution  are  important  enough  to  be  repeated  here,  not  only  to  see  how  a 
solution  is  effected,  but  also  to  gain  insight  into  the  nature  of  our 
model  as  a  servomechanism. 

Because  of  the  recursive  nature  of  equation  (2-l)  and  the  fact  that 
both  expressions  in  equation  (3-l)  are  recognizable  as  the  convolution  of 


5 


sequences,  it  is  clear  that  our  conditions  may  be  converted  by  means  of 
z-transforms  into  equivalent  and  more  algebraically  tractable  form.  To 


this  end,  if  V  represents  any  sequence  ,Vg,  •  •  • }  ve  vill  denote 

00  k 

the  corresponding  z-transform,  Y  V.  z  by  V(z).  Recalling  that  the 


k=0 


z-transform  of  the  sequence  (\+r)^_o  is>  (for  ®ny  integer  r)  z”rV(z) 
and  that  the  z-transform  of  the  convolution  of  two  sequences  is  given 
by  the  product  of  their  transforms,  ve  can  now  rewrite  equations  (2-l) 
and  (3-1)  respectively  as. 


(3-2) 


I(z)  =  zl(z)  +  z^+16(z)  -  X(z)  t  and 


(3-3) 


0(z)  =  A(z)  X(z)  +  B(z)  l(z) 


By  substituting  the  expression  for  0(z)  from  (3-3)  into  (3-2)  the 
inventory  transform  can  be  written. 


(3-4) 


1-z-z  ^(z) 


z^+^A(z)-l 

,  letting  S(z)  - - Srf -  >  ve  can  write,  more  simply, 


Or 


1-z-zT+1B(z) 


(3-5) 


I(z)  =  S(z)  X(z) 


thereby  expressing  the  inventory  at  any  time  t  as  a  convolution  of 
a  sequence  (S0,S1,Sg, . . . )  with  past  demands.  Because  of  (3-5)#  it  is 
quite  appropriate  to  call  S(  z )  a  transfer  function  and  to  view  inventory 
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as  the  output  of  a  "black  box"  the  input  of  vhich  is  demand  vhich  is 
operated  on  by  the  transfer  function.  Since  S(z)  depends  only  upon 
A(z)  and  B(z)  ve  can  define  the  latter  to  satisfy  our  criterion  of 
minimizing  expected  inventory  by  appealing  to  veil  known  results  from  the 
theory  of  servomechanisms. 


t 

Since  ve  can  write  1^  =  ^j^t-j  where 
t  t 

see  that  E(I^)  =  &  Sj  E<w  ■  £  Sj  n(t-j) 


Sj  is  non-random, 
for  t  =  1,2,3,  •  •  * 


ve 


Letting  E^.  =  E(l^)  denote  the  expected  inventory  at  time  t,  ve  see  that 

U(z)  =  X(z)  m(z)  in  terms  of  the  same  transfer  function  S(z).  Now, 

since  m(z)  is  deterministic,  the  requirement  lim  E.  =  0  is  equivalent 

t  -*  «> 

to  requiring  the  poles  of  s(z)  to  be  outside  the  unit  circle.  Since 

A(z)  only  occurs  in  the  numerator  of  S(z),  the  latter  requirement  depends 

only  upon  the  choice  of  B(z).  As  Vassian  suggests,  if  ve  let 

B(z)  =  -(l-zJ/Cl-z1*1),  then  S(z)  will  have  no  finite  poles.  Indeed, 

with  this  choice  of  B(z),  it  is  easy  to  see  that  S(z)  can  then  be 

/  T+l\  T 

written  as  (z^+1A(z)-l)  =  (z^+1A(z)-l)  Y  zk  vhich  becomes 

•  1-2  ra 


infinite  only  as  z  does.  We  have  thus  been  able  to  determine  one  of  our 
unknown  sequences,  though  the  choice  is  by  no  means  unique. 

As  to  the  determination  of  A(z)  ve  next  observe  that  ve  may  now 

write  l(z)  =  z1*1  A(z)  ?  X(z)  -  — -  X(z).  But 

1  T+l  T  . 

~  ^  zK  is  the  z-transform  of  the  sequence  (Uk  )£=()>  where 

Uk  =  1  if  k  <  T;  otherwise,  =  0.  Letting  A  (z)  =  z™A(z)  — <- 

*  /,  T+l\ 

ve  have,  from  equations  (3-4)  and  (3-3),  l(z)  =  A  (z)  -  - *-  x(z) 

1"Z 


X(z), 


7 


and  e(z)  =  — -  -1-^;  I(z).  Finally,  ve  may  write 
1MZ  z  1-z 

T+l 

(3-5)  «<*)  “  z‘(T+1V(z)  -  !(*)  >  a»4 

T+l 

(3-6)  A#(z)  =  I(z)  +  X(z) 


Or,  in  equivalent  form,  if  t  >  T, 


(3-7) 


j5)  =  At+l-HP  “  Xt  and 


(3-8) 


^t+l+T  "  Xt+1+T  + 


T 

j5)  Xt+1HJr-3 


t+l-KT 

w  +  jii  xj 


o 


Now,  from  (3-8)  E(At+1+T)  -  Et+1+T 


t+l+T 

+  T.  m(j).  Since  we  want 
j=t+l 


lim  E^.  =  0,  this  requires  that  we  choose  A^+1+T,  in  such  a  way  that 
t  — ♦  00 


E(A. 


t+l+T 


t+l+T 


)  -  Y,  m(J)”*0  as  t  -+».  Let  us  call  A^+1+T  a  forecast 


J=t+1 


of  the  demand  from  period  t+l  through  period  t+l+T,  i.e.,  over  the 
following  T+l  periods . 

We  have  thus  reduced  the  problem  to  one  of  determining  a  suitable 
means  of  forecasting  future  demands.  "Suitable"  here  means  that  the  forecast 
mujt  be  at  least  asymptotically  accurate  in  the  sense  of  the  above  limiting 
relation.  We  may  now  express  our  order  rule  in  the  form 
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g.wmw«Mwim 


ngfWM»B*i!  r-wwwH 


ttm9)  9t  =  At+1+T  “  0t-J  ’  Zt 

# 

vhere  A^+1+T  is  any  function  of  the  process  Xj.  such  that 

#  t+l+T  * 

E(A+.14~)  -  L  m(j)-40  as  t  -* «.  Of  course,  once  Al  is  specified 

*  j-t+l 

ve  can,  at  least  in  principle,  determine  the  sequence  Af  from  the 
definition  of  A  (z).  This  is  not  required  for  the  order  rule,  however, 
which  now  states  that  at  time  t,  one  orders  an  amount  equal  to  the  fore¬ 
cast  of  demand  over  the  next  T+l  units  of  time  less  those  orders  placed 
at  times  t-1,  t-2, ...,t-T  (which  will  be  receipts  at  times 

t+l,  t+2, ...,t+T)  and,  of  course,  less  the  inventory  at  time  t. 

All  that  remains  for  a  given  problem  is  to  determine  a  forecaster,  or 
equivalently,  the  function  A  (z).  No  further  explicit  expressions  can  be 

derived  until  more  assumptions  are  made  and,  quite  obviously,  the  deter- 
* 

mination  of  A  will  not  be  unique. 

4.  The  Model  as  a  Servomechanism 

The  general  inventory  system  which  we  have  described,  first  by 
equations  (2-l)  and  (3-l)  and  then,  in  a  transformed  version,  by  means  of 
equations  (3-2)  and  (3-3),  can  be  considered  a  discrete -variable  servo¬ 
mechanism.  Indeed,  as  we  have  brought  out  in  the  preceding  section, 
particularly  with  equation  (3-5)>  the  inventory  may  be  thought  of  as  the 
output  of  a  system  whose  input  is  demand.  In  fact,  in  the  present  case, 
the  system  resembles  an  automatic  control  device  with  a  feedback  system 
which  allows  for  taking  past  history  into  account  at  each  decision-making 
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stage.  These  facts  are  borne  out  more  clearly  by  standard  servomechanism 
diagrams.  Moreover,  such  diagrams  often  make  it  possible  to  replace  one 
system  by  an  equivalent  one  for  purposes  of  further  analysis. 

A  fev  remarks  concerning  the  meaning  attached  to  various  graphical 
forms  in  systems  analysis  may  be  helpful.  First  of  all,  vhenever  two 
transforms  are  connected  in  series,  we  will  alvays  understand  the  output 
to  be  the  product  of  the  two  quantities.  Thus,  the  basic  product  relation¬ 
ship  of  our  system  expressed  by  equation  (3-5)>  i.e.,  I(z)  =  S(z)X(z) 
can  be  portrayed  graphically  as  in  Fig.  4.1. 


Figure  4.1:  Graph  of  l(z)  =  S(z)x(z) 

On  the  other  hand,  whenever  two  quantities  are  connected  in  parallel, 
the  result,  usually  depicted  as  a  circle  or  node  at  the  junction  of  two 
directed  paths,  is  taken  to  be  the  sum  of  the  two  quantitites  as  shown 
in  Figure  4.2. 


Figure  4.2:  Graph  of  H(z)  =  F(z)  +  G(z) 
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With  these  conventions  in  mind,  the  inventory  system  described 
by  equations  (3-2)  and  (3-3)  can  be  represented  as  in  Figure  4.3*  In  this 
form,  the 


feedback  property  of  the  order  rule  is  dramatically  displayed. 

In  order  to  obtain  a  graphical  view  of  the  role  played  by  the  fore- 
*/  \ 

caster  A  (z;  of  our  last  section,  it  is  necessary  to  replace  the  original 
system  represented  by  Figure  4.3  by  an  equivalent  one.  This  is  accomplished 
by  means  of  repeated  applications  of  the  two  basic  rules  for  displaying  sums 
and  products,  respectively,  as  parallel  and  series  connections.  The  results 
are  displayed  in  Figure  4.4  which  displays  an  equivalent  system  free  of 


Figure  4.4:  Inventory  System  -with  Forecast 

Many  other  equivalent  systems  may  be  formed,  depending  on  the  various 
aspects  of  the  system  that  need  to  be  analyzed.  Also,  a  change  from  one 
form  of  pictorial  representation  to  another  often  suggests  definitions 
of  quantities  that  may  add  to  a  deeper  understanding  of  the  system. 

Once  a  flow  graph  of  a  system  is  determined  as  in  Figure  4.3  it  is 
possible  to  obtain  the  system  transfer  function  S(z)  by  means  of  a  few 
conventional  definitions  in  a  rather  easy  manner.  Indeed,  the  method 
often  is  a  distinct  advantage  over  the  analytical  solution  which  we 
determined  in  the  last  section. 

First  of  all,  we  shall  refer  to  any  function  denoted  by  a  square 
between  input  and  output  as  a  branch  transfer  function.  We  have  already 
referred  to  circled  functions  as  nodes  and  a  path  is  always  indicated  by 
arrows.  A  loop  is  a  closed  path  in  the  system  and  one  or  more  loops  always 
occur  in  feedback  systems.  For  any  simple  path  between  input  and  output, 
the  path  transmission,  P(z)  is  the  product  of  all  branch  transfer 
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functions  in  the  path.  For  a  loop,  the  product  of  all  branch  transfer 
functions  in  the  loop  is  called  the  loop  product.  Finally,  the  system 
determinant  A(z)  is  defined  to  be  unity  minus  the  sum  of  all  loop 
products  in  the  system.  Each  path  vith  transmission  Pj(z)  in  turn  has 
a  path  determinant  Aj(  z)  defined  as  unity  plus  the  loop  product  of  all 
loops  in  the  system  which  have  no  node  in  that  path.  With  these  definitions, 
it  turns  out  that  the  transfer  function  S(z)  can  be  easily  expressed 
by  means  of  the  equation, 

£Pj(z)  ^(z) 

(^-i)  S(z)  =  J ^ - 


In  the  present  case,  using  Figure  4.3,  there  is  a  single  loop  with 


product 


z^^z)  =  _  zT41 


1-z 


1-z 


T+l 


and  hence  the  system  determinant  is  given 


by  A(z)  *  1  + 


T+l 


1-z™  1-z™ 


.  There  are  two  simple  paths  with 


respective  path  transmissions  and  path  determinants  given  by, 
l  T+l. /  \ 

?l(.)  ■  '  ^(z)  =  1  and  =  ~  >  *2^  s  1*  Hence  we 


have  S(z)  *=  - -  1-z  - ^  as  ve  Previously  found. 


5.  Special  Cases 

One  of  the  primary  purposes  of  finding  the  system  transfer  function 
S(z)  as  discussed  in  Sections  3  and  4  is  to  be  able  to  examine  the 
system  response  to  various  demand  inputs.  Now,  with  the  inventory  system 
expressed  in  the  form  l(z)  =  S(z)  X(z),  once  it  is  assumed  that  the 
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demand  is  random,  it  follows  that  the  inventory  level  is  also  random. 

When  this  is  the  case,  it  is  of  course  not  possible  to  exercise  deterministic 
control  over  the  inventory  level.  For  one  thing,  there  can  be  no  explicit 
expression  for  the  demand  input  function  X(z)  in  this  case.  However, 
we  saw  in  Section  3  that  the  expected  inventory  level  responds  to 
average  demand  in  a  precisely  analagous  fashion,  l.e.,  E(z)  =  S(z)m(z) 
is  a  system  with  exactly  the  same  transfer  function  and  a  deterministic 
input  given  by  m(z),  the  z-transform  of  the  mean  value  function  of  the 
process.  With  this  in  mind,  we  now  proceed  to  examine  the  response  of 
expected  inventory  to  average  demand  in  some  special  instances. 

Certainly  one  of  the  simplest  cases  we  can  examine  is  that  of  a 
constant  mean  value,  i.e.,  we  suppose  that  m(t)  s  a  where  a  >  0.  For 

a 

such  a  case  then  m(z)  =  .  What  is  the  response  of  the  system  to  such 

an  input?  In  order  to  answer  the  question  it  is  necessary  to  know  the 
transfer  function  explicitly  and  this  in  turn  depends  upon  the  particular 
method  of  forecasting  adopted. 

Because  of  the  basic  feedback  property  of  the  system  we  have  defined, 
every  bit  of  past  history  is  available  whenever  a  decision  on  the  order 
quantity  must  be  made.  It  would  therefore  seem  natural  that  we  would  want 
any  method  of  forecasting  we  adopt  to  be  efficient  in  the  sense  of  utilizing 
all  of  the  information  available  at  a  given  order  time.  One  such  method  is 
that  of  exponential  smoothing  developed  somewhat  extensively  by  Brown 
[1]  and  others.  With  this  method  of  estimation,  some  credit  is  always 
given  to  each  past  demand  with  the  least  credit  going  to  the  oldest  demand. 
Moreover,  as  time  progresses,  the  credit  allowed  for  early  demands  becomes 
negligible. 
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If  ve  apply  exponential  smoothing  to  the  present  case,  then,  at  any¬ 
time  t,  ve  must  forecast  the  demand  over  the  next  T+l  periods. 

Since  our  best  estimate  of  demand  at  any  future  period  Is  given  by  mean 
demand  vhlch  in  turn  is  the  constant  a,  the  smoothed  estimate  of  a 

multiplied  by  T+l  will  do.  For  the  demand  history  X^,X^,. .  .,Xf.*  the 

t 

exponentially  smoothed  estimate  of  a  is  given  by  a  =  a  Y  p  X.  . 

k%>  Z~K 


vhere  a  is  the  chosen  smoothing  constant  (typically  0.1  to  0.3)  and 

*  ,  \  A 

p  =  1  -  a.  We  define  our  forecast  as  =  a  and  °^serve  that 

E(At+l-KT)  “  <T+1)  E(a)  =  (T+l)  a  £  BkE(Xt_k)  =  a(T+l)(l-pt+1)  and 
t+1+T 

J]  m(j)  =  (T+l)  a  in  the  present  case  so  that 
j=t+l 


e(a; 


*  t+i+T  t+l  * 

t+1+T)  “  ^  ^  “(J)  =  -(T+1)P  -+0  as  t-»».  Consequently, 


satisfies  our  basic  requirement  of  asymptotic  accuracy  for  a  forecast  as 
vas  established  in  Section  3. 

It  is  instructive  to  observe  that  in  the  present  case,  defining 

At+1+T  =  (T+1^  ^  is  e<luivalent  t0  letting  A(z)  =  • 

z  (l-pz)(l-z  ) 


Then  the  explicit  transfer  function  is  determined  as 

S(z)  =  ^  ^  '  From  servo  theory,  the  steady  state  behavior  of 

the  system  E(z)  =  S(z)  m(z)  is  defined  as  lim  E,  and  may  be  computed 

t  — » 

in  the  transform  domain  as  lim  (l-z)E(z).  In  the  present  case, 

z  -*  1 

(l-z)E(z)  =  S(z)  and,  since  lim  S(z)  =  (T+l)  -  (T+l)  =  0,  it 

z  -4  1 

follows  that  lim  (l-z)E(z)  =  0,  verifying  the  fact  that  lim  E.  =  0. 
z  -♦  1  t  ->  oo 

In  such  acase,  ve  say  there  is  no  steady  state  error. 
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As  a  second  example,  let  us  consider  the  problem  mentioned  in 
Section  1  of  accounting  for  changes  in  demand  when  an  item  changes  from 
a  slow  mover  to  a  moderately  fast  mover  to  eventually  a  fast  mover.  A 
first  approximation  to  describing  this  situation  in  the  present  model  is 
to  assume  that  m(t)  =  bt  for  some  small  positive  constant  b.  For  such 
an  item,  mean  demand  would  be  negligible  at  first  and  would  begin  to 
become  serious  at  some  future  time  depending  on  b. 

Of  course  we  do  not  know  the  value  of  b,  but  then  our  inventory  system 
is  set  up  in  such  a  way  as  to  respond  to  changes  in  demand  as  we  go  along  in 
time.  Indeed,  we  appeal  to  exponential  smoothing  once  more  as  a  means  of 
estimating  b  at  each  point  in  time  utilizing  all  of  the  past  data.  Other 
techniques  such  as  least  squares  might  also  be  used  to  estimate  b.  One 
of  the  reasons  we  prefer  exponential  smoothing  to  least  squares  is  the  fact 
that  far  less  information  has  to  be  stored  in  order  to  compute  a  running 
estimate  as  shown  in  Brown  [1]. 

For  the  present  case.  Brown  [l]  and  Reilly  [3]  have  shown  that,  for  a 

demand  history  of  X^,Xg,  ...,X^.,  an  appropriate  estimator  for  b,  based 

a  a  rA  ^  , 

on  exponential  smoothing,  is  given  by  b  =  -  [X^-X^j  where 


A  t  ,  ^  t  j^/v 

X,  =  a  Y  f3KX  and  X,  =  a  Y  p  X,  ,  .  Here  again,  a  is  the  smoothing 
t  ^  t-k  t  ^  t-k 

constant  and  X^  is  often  referred  to  as  the  double  exponential  smoothing 

_  ,A  A 

operator.  Reilly  [ 3 J  has  shown  that  lim  E(P;  =  b  so  that  B  is 

t  -*  00 

asymptotically  unbiased. 

t+l+T  A 

Because  of  the  above  results,  we  define  A,  ,  )  P  k  as  our 

t  1  1  k=t+l 


forecast  at  time  t  of  demand  over  the  next  T+l  periods.  Now,  for  t 

*  t+i+T 

sufficiently  large,  E(A,  +1+m)  =  b  £  k  and,  of  course 

z  1  L  k=t+l 
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t+l+T 


k*^ET+l 


m(k) 


t+l-HT 


k*J+l 


k  so  that  E(A. 


t+l+T 


t+l-HT 

)  -  J*  m(k)  -* 0  as  t  -» » 
k=t+l 


as  required.  If  ve  use  this  procedure,  ve  can  be  sure  once  again  that 
our  system  E(z)  =  S(z)  m(z)  has  no  steady  state  error.  Indeed,  the 
requirement  that  ve  adopt  an  asymptotically  accurate  forecaster  vas 
imposed  for  precisely  this  reason. 


6.  Summary  and  Recommendations 

We  believe  that  the  theory  of  servomechanisms  provides  some  powerful 
tools  for  analyzing  dynamic  inventory  systems.  The  present  report  has 
been  devoted  to  making  a  case  for  that  point  by  treating  one  particular 
model  as  a  discrete -variable  servomechanism  with  automatic  feedback.  Servo 
techniques  were  helpful  in  finding  explicit  /solutions  by  converting 
recursion  relations  to  z-transforms  where  the  mathematics  was  more  tractable. 
Even  more  important,  the  graphical  representations  are  extremely  helpful  in 
obtaining  an  overview  of  the  system  as  a  whole,  one  which  is  otherwise  often 
difficult  to  realize. 

There  are  several  directions  in  which  recommendations  for  further 
research  in  this  area  can  be  made.  As  previously  mentioned,  it  would  be 
of  great  practical  interest  to  process  a  real  inventory  system  in  which 
the  present  model  applies.  Using  real  data,  numerical  answers  could  then 
be  analyzed  toward  supporting  or  suggesting  changes  in  the  model.  Also, 
the  examples  treated  in  this  report  are  relatively  simple  ones  -  perhaps 
overs implified  for  practical  use  -  in  which  case  further  effort  should  be 

devoted  toward  applying  the  same  sort  of  analysis  to  more  complicated 
inputs . 


IT 


As  regards  the  model  Itself,  changes  can  alvays  be  suggested  that 
vould  typically  complicate  the  analysis  but  approach  reality  more  closely. 
For  example,  we  have  treated  the  lag  in  delivery  as  deterministic  whereas 
It  la  typically  randotn  in  reality.  The  model  could  be  changed  to  a 
continuous  review  system  in  which  case  z -transforms  would  be  replaced  by 
Laplace  transforms.  Presumably,  an  analysis  similar  to  the  one  we  have 
outlined  would  apply  to  that  case. 

Another  point  to  be  made  is  that  the  solutions  we  have  arrived  at 
are  explicit  only  up  to  specifying  a  forecast  of  future  demands  at  any 
given  time.  Now  there  are  many  ways  in  which  such  a  forecast  might  be 
specified  even  in  the  isolated  examples  we  have  treated.  Exponential 
smoothing  was  the  technique  adopted  in  this  report  for  the  particular 
reasons  mentioned.  Nevertheless,  other  methods  of  forecasting  might  well 
compare  more  favorably  on  further  reflection  of  the  desired  output  in  the 
model. 

We  have  imposed  the  restriction  that  forecasters  be  asymptotically 
accurate.  This  may  be  too  severe  and,  in  any  case,  one  without  that 
property  may  very  well  do  better  in  the  early  stages  than  one  meeting  the 
requirements.  If  such  forecasters  were  to  be  admitted,  then  explicit 
expressions  would  have  to  be  derived  for  the  inventory  level  at  any  given 
time  and  the  various  techniques  for  inverting  z-transforms  would  very  likely 
play  an  important  role  in  such  an  analysis. 

As  a  final  remark  it  should  be  observed  that  there  is  a  natural 
extension  of  the  servo  techniques  discussed  in  this  report  to  vector¬ 
valued  functions.  This  suggests  that  a  more  complex  multi-echelon  model 
might  possibly  be  treated  from  an  analogous  point  of  view  as  given  here. 
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Since  this  is  a  preliminary  report,  all  of  the  above  remarks  roty  be 
taken  as  suggestions  for  further  research  in  this  area. 
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